Determination of intergalactic magnetic fields from 

gamma ray data 



Warren Essey a , Shin'ichiro Ando b , Alexander Kusenko a ' c 

a Department of Physics and Astronomy, University of California, Los Angeles, CA 

90095-1547, USA 

b California Institute of Technology, Mail Code 350-17, Pasadena, CA 91125, USA 
C IPMU, University of Tokyo, Kashiwa, Chiba 277-8568, Japan 



Abstract 

We report a measurement of intergalactic magnetic fields using combined 
data from Atmospheric Cherenkov Telescopes and Fermi Gamma- Ray Space 
Telescope, based on the spectral data alone. If blazars are assumed to pro- 
duce both gamma rays and cosmic rays, the observed spectra are not sen- 
sitive to the intrinsic spectrum of the source, because, for a distant blazar, 
secondary photons produced along the line of sight dominate the signal. In 
this case, we find 1 x 10~ 17 G < B < 3 x 1CT 14 G. If one excludes the 
cosmic-ray component, the 10 -17 G lower limit remains, but the upper limit 
depends on the spectral properties of the source. We present the allowed 
ranges for a variety of model parameters. 

Intergalactic magnetic fields (IGMFs) play an important role in many 
astrophysical processes and may offer a new window on the cosmology in 
the early universe, but the size and origin of these fields are still poorly 
understood jlj]. Until recently, only the upper limits of 10 ~ 9 G were inferred 
from the observational data [2|. One can measure IGMFs using gamma-ray 
observations, for example, using time delays j^j, or by searching extended 
halos around the point objects [3, Here we employ an independent and 
complementary approach to obtain new upper and lower limits on IGMFs 
based on the spectra of three blazars observed by HESS, for which there 
also exist upper limits from Fermi Gamma-Ray Space Telescope. We have 
calculated the spectra numerically, and we use the goodness of fit to the 
data as the means to infer the average values of IGMFs. Blazar spectra 
have been used to set lower bounds on the IGMF 0, 0, S, @|- However, in 
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each of Refs. [si, 0, 0], only one EBL model and a single power law injection 
spectrum were considered, and the cosmic ray contribution was not included 
in Refs. JR 7j jl jjf , while it can be the dominant contribution for distant 
sources [nj, II, lil . Our analysis is different from previous work in that we 
have scanned over a very wide range of model parameters, and we included 
the cosmic-ray component. 

AGN are known to be sources of very high energy (VHE) gamma rays; 
they are also believed to produce cosmic rays 13|, |l4|. There is growing ev- 
idence that both primary gamma rays emitted at the source and secondary 
gamma rays produce d by cosmic rays along the line of sight contribute to 



the observed signals [10|, [III, Il2||. We, therefore, consider two cases: (i) 
pure gamma-ray emission at the source, and (ii) mixed gamma-ray and cos- 
mic ray emission. In the case (i) the intrinsic source spectra determine the 
observed spectra. If, however, the cosmic-ray contribution (ii) is included, 
the observed gamma-ray spectra are remarkably independent of the intrinsic 
gamma-ray or cosmic-ray spectra at the source. The observed spectrum is 
determined almost entirely by the shape (but not the overall normalization) 
of extragalactic background light (EBL) spectrum in case (ii), and the over- 
all agreement of such a prediction with observed spectra of distant blazars 
can be considered evidence of cosmic ray contribution, and, therefore, of 



cosmic ray acceleration in AGN [10|, llll, ll2|. Nevertheless, in what follows, 



we will consider both possibilities (i) and (ii). In the case (i), most of the 
observed secondary gamma rays are produced within tens of megaparsecs 
from the source, and our measurement reflects IGMFs closer to the source. 
In the case (ii), the field strength we infer corresponds to equally weighted 
averaging along the line of sight. 

To probe IGMF we constructed a detailed Monte Carlo simulation that 
accurately tracks the effects of both electromagnetic showers and cosmic- 
ray propagation and interactions. Included in the simulation are both pion 
photo-production and proton pair production for cosmic rays and pair pro- 
duction and inverse Compton scattering for gamma rays and electrons. Any 
relevant secondary decays such as neutron and pion decays were also in- 
cluded. Deflections due to the IGMF are tracked with both momentum and 
position vectors stored at all points of the propagation. To construct ob- 
served spectra, the position and momentum vectors of gamma rays at the 
z = surface are used to calculate the image size. Individual gamma-rays 
are tracked to the instrument, and their arrival direction further scattered 
by the instrumental PSF. Any events with arrival directions remaining with 
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the PSF radius are then retained. (A study of the source morphology as a 
function of energy can be a good diagnostic technique, but in this paper we 
concentrate only on the spectral analysis.) For E 1 < 100 GeV the publicly 
available Fermi PSF was used for comparison and for > 100 GeV HESS 
PSF was used. 

Suppression of the low-energy spectra of point images depends on the 
mean field along the line of sight more than on the (poorly known) distri- 
bution of IGMFs. Indeed, one can use the results of Ref. {§] to compare the 
effects of uniform magnetic fields with those of more realistic simulations, 



which take into account the distribution of matter in the universe [151 ] . The 
spectrum for "MHD(x)" ("Model 3x") model, with the mean field of 100 fG 
(0.1 fG), is close to that of the uniform 100 fG (0.1 fG) field. One can also 
understand this from general arguments 

0,1 IULIi2|- Thus, we modeled the 
IGMF as a characteristic strength organized into cubes of a characteristic 
correlation length each with a random direction. Of course, if the distribu- 
tion of matter along the line of sight is known, e.g., from Lyman-a forest 
in the direction to some blazars, one can further improve our inferences of 
IGMFs. 

From the available sources, we chose a subset based on three conditions: 
(a) the source must be observed at energies above 1 TeV, (b) the distance to 
the source must be relatively large, and (c) the source must show no temporal 
variability above 1 TeV. The latter condition is an necessary but not sufficient 
condition that secondary gamma rays dominate the observed signal, since 
primary gamma rays from blazars are predicted to show temporal variability. 
Based on these criteria, we chose three sources: 1ES 0229+200 (z = 0.14), 
1ES1101-232 (z = O.I86L1 and 1ES0347-121 (z = 0.188), each of which has 
been observed by HESS 16lll7lll8|. Temporal variability has been measured 
above the 200 GeV energy threshold of the HESS array for 1ES1101 and 
1ES0347 
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We consider a range of viable EBL models. A lower limit on the EBL 



density is set by galaxy counts [2l|. The upper limit is less known. Al- 



though gamma-raydata have been used to exclude higher EBL models such 



as Stecker et al. [22], it was recently pointed out that the inclusion of sec- 
ondary gamma rays from VHE gamma rays and cosmic rays significantly 
improves the fits to the data for these "high EBL" models [12j. Therefore, 
we include the model of Stecker et al. [22| as the upper limit (for UV EBL). 



A number of other models fall between these two limits 23] . EBL models 



with a higher level of UV background radiation are difficult to reconcile with 
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Figure 1: Photon spectra calculated numerically and normalized to fit HESS data points 
for 1ES0229+200 (z = 0.14) 0, 1ES1101-232 (z = 0.186) [Hj], and 1ES0347-121 (z = 
0.188) [3]. The Fermi upper limits at lower energy (thick solid line) were derived from the 
data in Ref. @. The three model predictions are for primary photons and B = 10~ 18 G 
(upper, dash-dotted line), B = 10~ 15 G (middle, dashed line), and B = 10~ 13 G (lower, 
thin solid line). The PSF of Fermi is folded into the signal below 100 GeV and the PSF 
of HESS is folded into the signal above 100 GeV. Here we assumed gamma-ray spectrum 
at the source oc E~ 1,75 with a 100 TeV high-energy cutoff. 
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a limit based on gamma-ray burst data [24 . 

We shall first discuss the case of primary gamma rays only. Following 
most recent papers on IGMFs, we will assume that the blazar spectra follow 
a power law dN/dE ~ E~ r with Y > 1.5 jiif. However harder spectra are 
also possible j26j. Furthermore, it is possible that the intrinsic spectra do not 
conform to a single power law. Deviations from single power law may affect 
the limits on IGMFs derived under the assumption of a power-law spectrum. 
We consider the range V = 0.5 — 2, which encompasses a broad variety of 
models. We also allow the high-energy cutoff of the power law spectrum to 
vary in a wide range, 10 — 10 4 TeV. 

We will consider magnetic fields ranging from 10 -18 G to 10~ n G. The 
IGMF correlation lengths are also not known; we assume it is in the range 
0.01-1 Mpc and use 1 Mpc as a reference point in Fig. 1. 

We fit the calculated spectrum to the data obtained by HESS Atmo- 
spheric Cherenkov Telescope (ACT) at high energy, perform statistical anal- 
ysis to calculate 95% confidence intervals, and we demand that the low-energy 
prediction not overshoot the Fermi upper limit. The spectra obtained for 
1ES 0229+200 (z = 0.14), 1ES 1101-232 (z = 0.186), and 1ES 0347-121 
= 0.188), are shown in Fig. 1 for IGMFs in the range 10~ 18 - 10~ 13 G, 
with a correlation length A c = 1 Mpc. The results for a source with intrin- 
sic spectrum T = 1.75, a 100 TeV high-energy cutoff and the "low" EBL 



model [21], shown in Fig. 1, illustrate the effects of IGMF on the spectrum. 
For larger IGMFs, more of the secondary gamma rays arrive outside the cone 
determined by the angular resolution of the instrument, with lower energy 
gamma rays arriving at larger angles. This causes an energy dependent drop 
in the spectrum. 

As one can see from comparing 10~ 15 G and 10~ 13 G curves in Fig. 1, 
the wash-out of the signal due to magnetic deflections causes a break in 
the spectrum at some energy, which is a function of the magnetic field. The 
following semi- quantitative analytical estimate may be helpful. The magnetic 



deflections depend on the energy [27 



e dcfl ~ 0.1°(1 + ^)- 2 (r e /5)- 1 ( J E; 7 /0.1 TeV)-\B /10~ 1B G), 

for A c 3> D c , where tq is the optical depth of the primary gamma ray and D e 
the electron cooling distance. For the energies considered here, tq oc K7 1 ' 2 , 
and for inverse- Compton (IC) cooling E^ oc E 2 , where E l0 is the energy of 
the primary photon, and n 7 oc E^ 1 / 2 , where n 7 is the number of photons 
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produced from IC cooling. Thus we expect G>d e fl oc B/E^ . Let us consider 
the flux of gamma rays that contribute to the point image (determined by 
a fixed instrumental resolution): 6d e fl < Qpsf implies E l0 > const x B 1 ^. 
Thus for an intrinsic power-law spectrum with index V we expect the point 
flux 

F oc — -7TT / — — ~ EZ 1/2 / Ej 



E] 12 JbU> dE l0 ' ~ ^ V/3 ^° ~ E ^ B (T-xy 3 ■ 

This relation implies the approximate scaling B 2 ( r ~^E 3 = const for the 
position of the break in the spectrum. By comparing the 10~ 15 G and 10~ 13 G 
curves in Fig. [1] we see that this approximation works reasonably well. Of 
course, this discussion is admittedly simplistic: the relative contributions of 
primary and secondary gamma rays, and some other effects were neglected. 
In deriving the actual limits we use the results of detailed numerical Monte 
Carlo for propagation, shower development, and the effects of magnetic fields. 

For angular resolution of about 0.1° we expect magnetic fields larger than 
about 1 fG to wash out the signal observed by the ACTs, while fields lower 
than a fG would boost the signal in the Fermi energy range. This allows 
us to use the combination of ACT and Fermi data to set both upper and 
lower limits on the IGMF. A lower limit can be found by requiring that the 
IGMF is large enough to wash out enough of the GeV signal so that the 
observed spectrum is below the Fermi 95% CL upper limits for the source. 
An upper limit can be found by requiring that the IGMF be small enough to 
avoid washing out the TeV signal. More specifically, we combined the ACT 
data for these three blazars and used 25 data points shown in Fig. 1. Then 
for each model with given B and T, we have three independent parameters, 
i.e., normalizations (or luminosities) of the spectra of three blazars. The 
number of degrees of freedom of the fit is 22. If the value of Xmin f° r each 
model specified with fixed B and T is larger than 33.9, then such a model 
is acceptable at the level of only 5% probability. This way, we set excluded 
regions (or upper limits on B) at 95% CL. 

Our results for pure gamma-ray sources are summarized in Fig. 2. The 
results depend on the choice of EBL model and high energy cutoff. We show 
both the "low" EBL model, based on lower limits from galaxy counts, and 
the "high" EBL model [22|, as well as high-energy cutoff values of 20 TeV 
and 100 TeV. 



Let us now include the the cosmic ray component [10j, [ll|, [12J . Now the 



observed gamma-ray spectrum does not depend on the intrinsic spectrum of 
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gamma rays, and it is also practically independent of the spectral properties 
of cosmic rays, such as the high energy cutoff and the cosmic ray spectral 
index (under some very mild assumptions) [l(J [HI 12]. Therefore, instead 



of the exclusion regions shown in Fig. 2, we find the following bounds (for 
A c = 1 Mpc): 

3 x 1(T 16 G < B < 3 x 1(T 14 G (High EBL) 

1 x 1(T 17 G < B < 8 x 1(T 16 G (Low EBL) (1) 

Of course it is possible, and likely, that the observed spectrum is a combina- 
tion of secondary particles from gamma rays and cosmic rays emitted by the 
source. In this case the most conservative allowed range for IGMF would be 
the union of the gamma-ray and cosmic-ray allowed ranges. 

The limits in Fig. 2 imply T < 1.8 for the set of blazars we chose, which 
is in good agreement with mean spectral index of 1.9 reported by Fermi for 
nearby sources. 

The limits presented in Fig. 2 are subject to some caveats. First, we 
have assumed that our three sources had similar intrinsic spectra, which 
need not be the case. Second, as discussed above, the likely contribution of 



protons along the line of sight [10|, llll, ll2| can reduce the limits to the model- 



independent form shown in Eq. [TJ Finally, time variability of the primary 
gamma ray sources can affect our conclusions [9j. 

Following the insightful discussion of Dermer et al. [9], it can be shown 
that the typical time delays for secondary photons are 

where Aioo ~ 1 for the Stecker et al EBL model, B_i 5 is the IGMF in units 
of 10~ 15 G. If the source is active for a timescale that is shorter than these 
delays, the observed signal is suppressed by a factor (t act i ve /At) where t a ctivc 
is the time that the source was active for. This suppression mimics the effects 
of a stronger magnetic field, and, therefore, the predicted limit on the IGMF 
would be suppressed by the a factor 



iartivs \ /l0 15 G\/ E, 



K 




1 



2 x 10 6 yrs J \ B limit J VlO GeV/ 
where -Bii m ;t is either the upper or lower limit shown in Fig. 2 



V (3) 
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Figure 2: Limits on the magnitude of the IGMF and intrinsic spectral index based on the 
combined fit for spectra of 1ES 0229+200, 1ES 1101-232, and 1ES 0347-121, assuming 
that only gamma rays are produced at the source (no cosmic ray contribution). Shaded 
regions are excluded at 95% CL from either Fermi or ACT data. The limits are shown for 
"low" and "high" EBL models and for two values of the high-energy cutoff: 100 TeV and 
20 TeV. IGMF correlation length is assumed to be 1 Mpc. 
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The observations of 1ES 0229+100 by HESS and VERITAS 
were separated by roughly 3 years and show no evidence of significant vari- 
ability. Variability may be present during the times when the source was 
not observed, or at a level below the instrumental sensitivity We assume 
^active ~ 3 yrs as a lower limit on the time the source was active. The most 
restrictive energy for the lower limits is near E~ ~ 10 GeV and most of the 
lower limits fall in the range 10~ 15 — 10~ 17 G which means the most conser- 
vative lower limits would be suppressed by a factor of 10 — 10 3 . For upper 
limits, the most restrictive energy is at E 1 ~ 1 TeV and since most upper 
limits fell into the range 10~ 13 — 10~ 15 G we expect the upper limits to also 
be lowered by roughly 10 — 10 3 . Thus, the overall effect of including a short 
activity time is to lower the 95% confidence intervals shown in Fig. 2 and 
our robust 10~ 17 G lower limit would be lowered to 10 -18 G in the most 
conservative case. 

It is important to distinguish between the time scale of the source ac- 
tivity and short-term variability. A variability of the source would not lead 
to a suppression discussed above because signals of variable activity would 
undergo significant delays in the IGMF and would be observed as a constant 
source with the time averaged luminosity of the flares. Thus, if the source 
is active on a timescale of 10 6 yrs, the limits in Fig. 2 are valid, even in the 
case of variability on shorter time scales. 

The effect the source activity-dormancy time scale is much weaker in the 
case of the cosmic ray limits, because the delays in the stronger magnetic 
fields near the source can wash out variability on even longer time scales. 
The magnetic fields within the host galaxies containing source are likely to 
be of the order of 1 fiG which can lead to significant delays. Fig. [3] shows the 
delays for a 100 kpc wide source of 10 10 GeV protons with B = 1 /iG and 
Z c = 0.1 kpc. It can be seen that the delays are significantly longer than the 
observed variability of blazars and thus we expect any intrinsic variability to 
be washed out. Such a source appears as a constant, long-lived source, and 
any time delays incurred in the host galaxy are not measurable. Of course, 
there is a significant uncertainty in the magnetic field distribution around 
the source, in particular, in the direction of the jet, and the values assumed 
here are taken solely to illustrate the idea. 

The effect of the source variability would be to damp the observed power 
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10 50 100 500 1000 5000 lxlO 4 

Arrival Time (yrs) 

Figure 3: Arrival time probability distribution in arbitrary units for primary cosmic rays 
including only delays at the source. A cosmic ray source of 10 10 GeV protons was used in 
a source 100 kpc wide with B = 10~ 6 G and l c = 0.1 kpc. 



of the source by a factor of 

f damp r " 1 -^active ~7 > (^0 

''delay 

where tdeiay is the typical proton delay at the source, t ac tive is the typical time 
the source is active or flaring and A^ active is the number of times the source is 
active in the time period ideiay This damping will not be a significant effect, 
especially since the typical deflections at the source will not be enough to 
affect the beaming factors. 

Using Eq. fl3]) and these typical time delays we see that the model inde- 
pendent lower limit of 10 -17 G would remain unchanged and the upper limit 
of 3 x 10 -14 G would be lowered by roughly one order of magnitude. Thus 
the limits we provided for cosmic rays are still the most conservative limits 
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even including all timing uncertainties. Thus, for the case of cosmic rays 
dominating the signal, we can report a model-independent 95% CL interval 
for the IGMF of 

1 x 1(T 17 G < B < 3 x 1(T 14 G. (5) 

Our results represent the first 95% confidence interval measurements of 
the IGMF based on the spectral data alone. The method used here to probe 
the IGMF will become more powerful in the future. Future observations 
of distant blazars will not only improve the limits on the IGMF, but also 
constrain AGN and EBL models. The measured values of IGMFs can be 
used to distinguish between the cosmological and astrophysical scenarios for 
the origin of all astrophysical magnetic fields. In addition, this information 
will lead to a better understanding of gamma rays and cosmic rays, as well 
as the properties of AGN and of universal photon backgrounds. 
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